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Abstract 

Causal "superluminal" effects have recently been observed and discussed in 
various contexts. The question arises whether such effects could be observed with 
extremely weak pulses, and what would prevent the observation of an "optical 
tachyon." Aharonov, Reznik, and Stern (ARS) [Phys. Rev. Lett. 81, 2190 (1998)] 
have argued that quantum noise will preclude the observation of a superluminal 
group velocity when the pulse consists of one or a few photons. In this paper we 
reconsider this question both in a general framework and in the specific example, 
suggested by Chiao, Kozhekin, and Kurizki [Phys. Rev. Lett. 77, 1254 (1996)], 
of off- resonant, short-pulse propagation in an optical amplifier. We derive in the 
case of the amplifier a signal-to-noise ratio that is consistent with the general ARS 
conclusions when we impose their criteria for distinguishing between superluminal 
propagation and propagation at the speed c. However, results consistent with the 
semiclassical arguments of CKK are obtained if weaker criteria are imposed, in 
which case the signal can exceed the noise without being "exponentially large." 
We show that the quantum fluctuations of the fleld considered by ARS are closely 
related to superfluorescence noise. More generally we consider the implications 
of unitarity for superluminal propagation and quantum noise and study, in addi- 
tion to the complete and truncated wavepackets considered by ARS, the residual 
wavepacket formed by their difference. This leads to the conclusion that the noise 
is mostly luminal and delayed with respect to the superluminal signal. In the limit 
of a very weak incident signal pulse, the superluminal signal will be dominated by 
the noise part, and the signal-to-noise ratio will therefore be very small. 
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1 Introduction 



Chiao et al. |T]],|0,|§] have shown that certain "superluminal" effects are possible 
without violation of standard notions of Einstein causality, i.e., without conveying 
information faster than the velocity c of light in vacuum. Such effects have been 
demonstrated experimentally in optical tunneling and in an electric circuit 

0- 

It has been suggested by Chiao, Kozhekin, and Kurizki (CKK) that an 
optical pulse can propagate superluminally in an amplifier whose relaxation times 
are long compared with the pulse duration. The dispersion relation they derive 
can be obtained directly, as follows, starting from the formula for the refractive 
index of a monatomic gas: 

n[u}) = l + 2^2^— 2 ^ 

m i j ^ji — ^ 

for n{uj) = 1, where iVj is the number density of atoms in state i and f{i,j) is the 
oscillator strength for absorption on the i ^ j transition of frequency Uji. Near a 
two-level resonance this becomes 

n{uj) = 1 + , (2) 

m uj^ — 00"^ 

where 1 and 2 designate the lower- and upper-energy levels, respectively, and 
t^o = ^21- Close to the transition resonance frequency Uo, 

nH = l + — ^ \ 1- 3 

mijJo — — 

when we include a dipole damping rate /3. The (real) refractive index near a 
resonance is then 

Introducing the inversion w = {N2 — Ni)/N, where is the number density of 
atoms, and assuming a field sufficiently far from resonance that {ujo — » 
we have 

ne^Nwf 1 

nR{uj} = 1 , (5) 



UJ UJ 

nR[uj)- = - 
c c 



1 - 



Tie^Nwf 1 
mujo uJo — ^ 



UJ 

c 



UJ. 



pW/iuJo 
UJo — UJ 



(6) 



, , 1. ^ uj uj^w/AuJo ^ 1 ujIw/4c 

k-ko= -[uj - uJo) = -[oj - uJo) , (7) 

C C UJq — UJ C UJr, — Ul 



and 



fi2 - Ken + -wujI 
4 ^ 
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where K = k — ko, = u — Uo, and the "plasma frequency" Up is defined by 

ujI = ATTNe^f/m = SuNSujJh , (9) 



with d the electric dipole transition moment. Equation (|^) is the dispersion relation 
obtained by CKK. 

We refer the reader to the CKK paper for a discussion of this dispersion relation. 
Here we simply note that implies the group velocity 



duj 



dk 



1 - 



ujIw/4: 



(10) 



so that, in the case of an amplifier (w > 0), a short off-resonant pulse can propagate 
with a group velocity Vg > c. 

Questions have been raised about the validity of the latter prediction at the 
one-photon level, which would correspond to what CKK call an "optical tachyon" 
[0. Aharonov, Reznik, and Stern (ARS) have presented general arguments, 
based on the unitary evolution of the state vector, that "strongly questions the 
possibility that these systems may have tachyonlike quasiparticle excitations made 
up of a small number of photons." They also consider a particular model as an 
analog of the CKK system. 

In this paper we address the question of superluminal propagation at the one- 
or few-photon level, and in particular the role played by quantum noise in the 
propagation of such extremely weak pulses. We begin in the following section with 
some physical considerations about the observability of superluminal propagation, 
and we briefly compare the ARS and CKK models. In Section 3 we formulate 
the Heisenberg equations of motion for the propagation of a short optical pulse 
in an inverted medium, and briefly review some relevant results from the theory 
of superfluorescence (SF). In Section 4 we derive a signal-to-noise ratio for the 
case where an incident, Gaussian signal pulse made up of q photons is very short 
compared with the radiative lifetime and has a central frequency far removed from 
the resonance frequency of the medium. If we impose the ARS criterion for dis- 
tinguishing between superluminal propagation and propagation at the speed of 
light, we flnd, consistent with their conclusions, that the signal must be "exponen- 
tially large" in order to distinguish it from quantum noise. If the ARS criterion 
is replaced by a much weaker one, however, the signal-to-noise ratio can exceed 
unity even for a one-photon signal pulse, as suggested by CKK. We relate the am- 
plifled quantum fleld fluctuations of ARS to quantum fluctuations of the atomic 
dipoles in the case of the optical amplifler. In Section 5, following the ideas of 
ARS, we present some general considerations based on the premises of unitarity 
and superluminal propagation. ARS show that, when the group velocity exceeds 
the speed of light, the superluminal signal is reconstructed from a truncated initial 
wavepacket, and that this truncated wavepacket has unstable modes. We show 
that the truncated wavepacket introduced by ARS propagates with both luminal 
and superluminal parts, and that, while the superluminal part is the reconstructed 
signal, the luminal part has the exponentially growing parts corresponding to the 
unstable modes. In addition, we study the residual wavepacket formed by the 
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difference of the complete and truncated wavepackets. We show that contribu- 
tions from the truncated and residual wavepackets cancel in the luminal region, 
but that, unlike the signal, the noise does not cancel, leading to the conclusion 
that the quantum noise is mostly luminal rather than superluminal. In the limit 
of a very weak incident signal pulse the signal-to-noise ratio will be very small, 
consistent with the conclusions reached by ARS. 

It may be worth recalling that a primary reason for rejecting the possibility 
of superluminal transmission of information is the requirement that causality be 
maintained when Lorentz transformations are made: superluminal transmission of 
information would allow an event A causing an event B in one reference frame to 
occur after event B in a different frame. Considerations of superluminal propaga- 
tion therefore often raise questions relating to Lorentz invariance. When and how 
should one include relativistic effects in order to ensure that physically meaningful 
results are obtained? 

As in all previous treatments of pulse propagation in an inverted medium that 
we know of, we choose the reference frame in which the atoms are at rest. The 
Lorentz invariance of the fundamental, fully relativistic theory implies, of course, 
that our conclusions do not depend on this specific choice of a reference frame. 
Working in this frame, we treat the response of the atoms to the field in the ap- 
proximation of nonrelativistic quantum mechanics. The electromagnetic field in 
this frame is also treated approximately, namely in the slowly varying envelope ap- 
proximation that is used practically universally in the theory of resonant atom-field 
interactions. A different choice of reference frame would require us to start with 
the fully Lorentz-invariant equations and then make the slowly-varying-envelope 
and other approximations as appropriate. These approximations are known to be 
very accurate unless, for instance, the light pulse is extremely short, and to the 
extent that they are valid our results and conclusions are Lorentz-invariant. 

2 Preliminary Considerations 

The quantum noise limitations to superluminal propagation discussed by ARS were 
associated physically with spontaneous emission in the case of an optical amplifier, 
and could invalidate the CKK results in two ways. First, CKK assume that the 
atoms stay in their excited states as the pulse propagates through the amplifier. 
Radiative decay of the excited state will modify their "tachyonic dispersion re- 
lation" and, if the decay is rapid enough, can lead to a subluminal rather than 
superluminal group velocity, since w in equation (|10D can become negative. This 
can be avoided by using a sufficiently short pulse. Second, spontaneously emitted 
radiation might interfere with the measurement of the superluminal group velocity 
by introducing substantial noise. It is this possibility that is addressed by ARS. 

Although the ARS arguments are certainly compelling, they are based in part 
on an analog of an optical amplifier rather than a theory involving the interaction 
of the electromagnetic field with an atomic medium. In particular, theirs is a model 
of a single quantum field rather than coupled atomic and electromagnetic quantum 
fields. The dispersion relation associated with this model, and the criteria assumed 
by ARS for the observability of superluminal propagation, lead to the conclusion. 
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by analogy to an optical amplifier, that spontaneous emission noise cannot be 
avoided no matter how short the pulse or the transit time through the amplifier. 
Specifically, the unstable modes appearing in their model - which "are analogous 
to spontaneous emission in the optical model of an inverted medium of two-level 
systems" |^ - will preclude the observation of superluminal group velocity when 
the pulse is made up of a small number of photons; the quantum noise will be 
larger than the signal. In this section we present some physical considerations, 
motivated by the CKK and ARS analyses, for the observability of superluminal 
group velocity. 

Following their equation (11), ARS state two necessary conditions for the ob- 
servability of superluminal propagation [c = 1 in their units] : 

(1) VgT » l/6k, where Vg is the group velocity, T is the time at which the 
wavepacket is observed, and 5k is the spectral width of their initial pulse. 

(2) {vg - 1)T » l/6k . 

The first condition ensures that "the point of observation [is] far outside the initial 
spread of the wavepacket." The second allows us to "distinguish between superlu- 
minal propagation and propagation at the speed of light." 
In the ARS model, where the field (p satisfies 

_|__|_,,„.^ = 0. (11) 
= , , (12) 



the group velocity is 



-'9 

m 



^2 no 2 



where ko is the central value of the spatial frequency k for the initial pulse. For 
m < ko we can approximate f ^ by 1 -|- m^/2/c^, so that condition 2 [and also 
condition 1] is satisfied if 

m^T » kl/5k » ko . (13) 

ko >> 1/T - the condition that the observation time should be much larger than 
the central frequency of the pulse - then implies 

mT»l . (14) 

Since for mT >> 1 the amplified quantum noise grows exponentially [see Section 
3], ARS conclude that the "signal amplitude should be exponentially large" in 
order to distinguish it from noise. Thus, according to ARS, the observability 
of superluminality for an input pulse consisting of only a few photons would be 
clouded by spontaneous emission noise. 

Consider now the implications of conditions 1 and 2 for the actual system of 
interest, namely a very short optical pulse in an inverted medium. Can we satisfy 
these conditions for observation times short compared with the radiative lifetime? 
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For a short optical pulse of central frequency uj propagating in an inverted 
medium {w = 1) with resonance frequency Uo, the refractive index is [equation 



2'KN(P/h 



2 



n{u) ^ 1 + '— = 1 - (15) 

for Up/{4:Uo) « \u!o~uj\ = |A|. We are assuming that |A| is large compared with 
the absorption width, which in our case is the radiative decay rate. Equation (|15D 
implies 

!!l=( -^[coniio)]] = ]—— (16) 



and 

c 1 - CU2/4A2 4A2 c ■ ^ ' 

Then conditions 1 and 2 of ARS become, respectively. 



>> 



1 - cu2/4A2 c6k ' ' 

with Tp the pulse duration. Both conditions can be satisfied if, for instance, 
T » Tp and uTp/AA"^ is not too small. To avoid spontaneous emission during 
the observation time T, take T « tj^ad, where trad is the radiative lifetime of a 
single inverted atom. Then the ARS conditions require that 

Trad » T»Tp. (20) 

As noted by CKK, there is another aspect of an inverted atomic medium that 
must be addressed, namely superfluorescence (SF). SF is a collective phenomenon 
of the sample as a whole. We shall denote by Nt, S, and L the number of atoms, 
the cross-sectional area, and the length of the sample, respectively, so that the 
density of atoms is given by = Nt/SL. If coUisional and other dephasing 
mechanisms are sufficiently weak, an inverted medium of Nt atoms can emit SF 
radiation at the rate 

tr = trad/Nt , (21) 

i.e., the radiative decay time can in effect be smaller by a factor of than the 
single-atom radiative lifetime trad assumed in the discussion thus far. The peak 
of the SF pulse occurs at a time 



^ ln(27rA"T) 



2 



(22) 



following the excitation of the atoms. It would appear then that the quantum 
noise associated with SF will be small if 

Tp, L/c <tr<td . (23) 
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We note for later purposes that 



^p = 7- = -wtSc=—-, (24) 

^ n tradSL trL 



where we have used equation ( 131| ) of Appendix A for the single-atom radiative 
lifetime trad- 

This brief summary lends support to the CKK suggestion, but obviously a more 
quantitative analysis is called for. To this end we now formulate, in the Heisenberg 
picture, the quantum theory of pulse propagation in an amplifier. 



3 Formalism for Pulse Propagation 

We begin with the Hamiltonian for Nt two-level atoms (TLAs) interacting with 
the quantized electromagnetic field via electric dipole transitions: 

H = -huo CTzj - dYa^jE{zj) + ^ hujkd\dk , (25) 

j=l j=l k 

where uJq and d have the same meaning as before and Zj is the z-coordinate of atom 
j. The carets ( " ) are used to denote operators. We consider a one-dimensional 
model in which the atoms occupy the region from z = to z = L and the field is 
a superposition of plane waves propagating in the z direction. The electric field 
operator is given by E{z) = E^^\z) + E^~\z), where 



k 



^k = u,/c) (26) 



and E^~\z) = E^'^^zY . Si, where S, as before, is a cross-sectional area and 
i a length, is the quantization volume. For simplicity we consider only a single 
field polarization, namely linear polarization along the direction of the transition 
dipole moment of the TLAs. dk and al are the photon annihilation and creation 
operators, respectively, for mode k, and the a's are the Pauli two-state operators 
in the standard notation. 

We will work in the Heisenberg picture, in which the time-dependent electric 
field operator satisfies 

f I ^ 47r d^P And ^ O^a^^ , 47r , 9^ . , , 

(27) 

where in the last step we have made the continuum approximation for the polar- 
ization density P, assuming a uniform atomic density A^. We now write 

i W {z, t) = F{z, t)e-^'^(*-^/^) (28) 
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and assume F{z,t) is slowly varying in z and t compared with exp[—iu{t — z/c)]. 
In this approximation 

It will be convenient to use the atomic lowering and raising operators a = ^{ax — 
icTy) and = \{crx + icFy)^ respectively, such that [a, a'^] = — a^, and to write 

a{z,t) = s{z,t)e-'''^'-'/^\ (30) 

where the operator s{z,t) is assumed to be slowly varying in the same sense as 
F{z,t). Then, in the rotating-wave approximation (RWA), we can replace ( ^91) 
with 

f + if = (2.,m^)s . (31) 

OZ C Ot C 

where on the right-hand side we have approximated u by Uo- This equation and 
the TLA Heisenberg equations 

f^=-^{A~^f3)s-Ja.F, (32) 

^ = -2/5(l + a.)-|^(FtS-StF), (33) 

derived in Appendix A form a closed set of operator equations. They provide the 
basis for a quantum theory of propagation in either amplifying or absorbing media. 

In the semiclasical approximation in which the atom and field operators are 
replaced by their expectation values, equations (pl|)- (|3^) reduce to well known 
Maxwell-Bloch equations. Otherwise, different limits can apply: 

• The limit of /? ^ 0, A = and cr^ ^ 1 considered below gives equations 
(|35|)-(|37D implying superfluorescence when the initial state of the field is the 
vacuum. 

• The limit of ^ cuq gives the ARS field equation, as discussed below. 

• Finally, in Section 4 the CKK case of large detuning, 6"^ ^ 1, and the initial 
state of a very short incoming pulse, is studied. 

If the field central frequency u is assumed to match exactly the atomic res- 
onance frequency Uq, so that A = 0, and if we restrict ourselves to times short 
compared with the single-atom radiative lifetime [trad = (2/3)"^] and assume that 
the atoms remain with probability = 1 in their excited states over times of interest, 
we can ignore ( ^3|) and replace crz{z,t) by 1 and equation (|32|) by 
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In terms of the independent variables C, = t — z/ c and t] 



implying 



^ = ( 2mNd'^] s , (36) 

OT] \ C J 



dr]dC \4c/ dr]dC \ 4c 



Equations (p5|)-(^) have been used in studies of the buildup of superfluorescent 
radiation [QT^It will be useful for the discussion in Section 4 to briefly rederive 
here one of the most important results of those studies. 
Equation (pT]) has the formal solution 



F{z,t) = F,{z,t) + (2mNd'^] r dz's{z',t-^^)e{t ^ ^' ' 

\ c J Jo c 



c 



= Fo{z, t) + \ 2mNd^j dz's{z -z',t- z'/c)e{t - z'/c) , (38) 

where we have chosen the retarded Green function over the advanced Green func- 
tion in order to ensure causality. Here 9 is the unit step function and Fo{z,t) is a 
solution of the homogeneous equation. We are interested here in the expectation 
value {F'^{L,t)F{L,t)) at the end {z = L) of the medium. For SF the expectation 
value is taken over the vacuum state of the field, in which case the first term on the 
right-hand side of (|38D does not contribute to normally ordered expectation values. 
We may therefore ignore this term for practical purposes. Defining y = we 
find from (|37|) that s satisfies the differential equation for Io{y), the modified Bessel 
function of order zero |jlO|- The solution of interest for F{L,t) is then |^ 

F{L,t) = (2TTiNd'^^ £ dz's{L-z\0)Io [up^ {z' / c){t - z' / c)^ e{t-z'/c) . (39) 

In order to calculate {F\L,t)F{L,t)) we require (s''"(2;', 0)s(z, 0)), which we 
evaluate in Appendix B. We obtain 

{F\L,t)F{L,t)) = (^2nd^y^£dxeit-x/c)I^o (^P\/(^/c)(t - x/c)) . (40) 
For times large enough that Iq may be replaced by its asymptotic form, 

(Ft(L,t)F(L,t))~^^eV^. (41) 
Stt Set 

Equating the intensity expectation value (c/27r)(F^(L, t)F(L, t)) to the maximum 
expected SF intensity NtTi-uJo/ St^, we arrive at the expression (E2[) for the time 
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at which the SF pulse reaches its peak intensity. In the short-time hmit, on the 
other hand, 

/ , . \2 AT 

{F\L,t)FiL,t)) [2nd^j -ct , (42) 
a result we will return to in Section 4. 



Approximation Leading to ARS Field Equation 

Our considerations thus far assume that the field central frequency lies in the 
vicinity of the atomic resonance in the sense that the detuning A is small in 
magnitude compared with u and Uo- Let us now suppose instead that the field 
frequency u is very large compared with Uo- In this case we must work with the 
atomic operators ax, CFy instead of the slowly varying s. /^From equations 
and ( |121| ) of Appendix A we have 

a, + uola, = T-^a,E ^ -^E 43 

n h 

in the approximation = 1. The assumption uj » implies 

- (44) 



0", 



so that, from equation ( P7| ) 




(45) 



This is identical to the equation of motion for the quantum field in the ARS model 
when we equate to their m^. From this perspective the ARS equation of motion 
describes the interaction of the electromagnetic field with unbound electrons 
{u » Uo) per unit volume. However, the usual plasma dispersion formula = 
1 — uijio^ for the refractive index n is replaced in this case by 



n 



,2 1 , , .2 /, ,2 



ujIIuj' . (46) 

This is a consequence of the assumption = 1; had we assumed = — 1 we 
would have obtained the familiar plasma dispersion formula. 

To describe the growth of the quantum noise with time in this model, we write 
( ^5]) in the form 

where ti = t — z / c, T2 = t + z / c. In terms of the independent variable y = my/r^, 
equation (^71) has solutions that are linear combinations of the zero-order modified 
Bessel functions Io{y),KQ{y). For large t, the vacuum expectation value 

„2mt 

SO that the quantum noise grows exponentially in time from the initial fiuctuations 
of the vacuum field, the fiuctuations present before the medium in the ARS model 
is "inverted." 
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4 Signal and Noise 



We wish to determine to what extent the observation of the superluminal group 
velocity considered by CKK will be affected by quantum noise. The system of 
interest is described by the Heisenberg equations of motion (|3l|) and (|3^ . We 
approximate Oz by 1, assuming that pulse durations and transit times Ljc are 
sufficiently small that de-excitation of the initially inverted atoms by radiation (or 
any other decay process) is negligible. The situation here is different from that 
describing the onset of SF in that (a) the detuning A is not zero but is instead large 
(Section 2), and (b) the initial state of the field is not the vacuum but corresponds 
to a short pulse of radiation from some external source. 
The equation for 5(2;, t) in the present model is 

| = -.(A-z/3)S-^F, (49) 

or 

S(Z, t) = S{Z, t„)e-(A-/3){t-to) f ^t'P^^^ ^>yiA-im'-t) . (59) 

n Jto 

to is some initial time, before any pulse is injected into the medium. We take 
F{z,to) = 0, although of course what this really means is that there is no non- 
vanishing field or intensity in the medium at to, so that for practical purposes 
(normally ordered expectation values) we can in effect ignore the operator F{z, to) 
in the equation for s{z,t). 

The pulse is assumed to have a central frequency u and to have no significant 
frequency components near the resonance frequency uJo'- l^kp > 1- We assume 
that |A|rp is large enough that we can approximate (^) by integrating by parts 
and retaining only the leading terms: 

- *(M.)e--'-«<-> - |A±|^(M) - . (51, 

As will be clear from the analysis that follows, this approximation implies the 
undistorted propagation of the incident pulse at the group velocity Vg, as assumed 
by CKK. 

/^From (H), 



l^+i^ ^ (2mNd^) S(z,ge-(^-/^)(-°)+|F+.[n(o;)-l]5^F+fl - ^ 
az c ot \ c J 2 c \c VgJ ot 

(52) 

where 

ATI Nd'' UJo /3 
' = he A^ + f3^ ^'^^ 

is the gain coefficient for propagation of a field with frequency a; in the inverted 
medium. We have used equation (|I5|) for the refractive index n{uj) and {\rf\j for 

Vg/c- 1. Writing F{z,t) = F'(z, t)e^["('^)~i]^^/" and s{z,to) = s'{z,to)e'^''^'^^-'^^^'^'' 
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yields an equation in terms of the primed variables in which the term i[n{u!) — 
l]{u}/c)z associated with phase velocity is eliminated. Then, ignoring for practical 
purposes the difference between the primed and unprimed variables, we have 



and therefore 

F{z, t) = F(0, t - z/vg)e''/^ + (^2mNd'^^ dz's{z\ Qe''^''''^'^ 

^ g-*{A-i/3)[t-t„-(.-.')/.,]^(^ _ _ _ z') /Vg) 

= F,{0,t-z/vg)e^'/^ + Fn{z,t) , (55) 
where the subscripts s and n denoted signal and noise, respectively. Here 

Fn{z,t) = (2mNd'^^ j[^'rf2's(z',to)e3(^-"')/'e-'(^-''^)[*-*°-("-^')/'''']0(t-to-(2-2')K) 

(56) 

is a quantum noise field associated with the quantum fluctuations of the atomic 
dipoles. 

To appreciate the significance of g as defined by equation (|53|) , consider the gain 
coefficient qr for a radiatively broadened transition of frequency Uo and radiative 
decay rate I/trad = 2/5. For light of frequency u = Uo — A, 

_ NS 2(3 _ A-KNd^uJo (3 

tradA^ + (3^ ~ he A2 + /52 ^ ^ 

if we assume that all the N atoms per unit volume are in the upper state of 
the amplifying transition. Thus qr = g, i.e., g is just the gain coefficient for 
amplification by stimulated emission. We note also that, from equation ([T7|), 



g = 2/3 ---] {5t 



^9, 



in the case under consideration where the amplifying transition is radiatively broad- 
ened and the detuning is large compared with the gain bandwidth. 

The operator s{z,to) has the expectation-value properties ( |140| ) and ( |141| ) of 

Appendix B. These properties imply {Fn{z,t)) = {F^{z,t)) = and 



{Fl{z,t)F^{z,t)) = (2nNd^) ^e'^^^*"*") f d/e^(— ') 

V C / I\T Jz-Vg{t-to) 

(59) 



c , 



\2 AT c 


' (.9^ai _ e-2/3t" 


/ 'S2^ 





where we have used the relations (|58D and = NSL and, to simplify the notation, 
we have taken to = 0. 
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Since the atom and field are initially uncorrelated, i.e., 

{F\0,t-z/vg)s,ito)) = {FliO,t-z/vg)){s,iQ) = 0, (60) 
we have, at the end of the amplifier, 

{F\L,t)FiL,t)) = {F}iO,t- L/vg)FsiO,t- L/vg))e^'^ + {FliL,t)F4L,t)) (61) 
and the signal-to-noise ratio 



SNR{L,t) 



{FjiO,t-L/vg)F,iO,t-L/vg))ea^ 
{Fl{L,t)F4L,t)) 
{F}{0,t-L/vg)F,{0,t-L/vg))e^^ 
{2TTdujo/cY{N/S){c/2(3) [e9L - e-2/3i/-s] 

(Ft(0, t - L/vg)Fs{0, t - L/t;g))e^^(^/^-^A'^)^ 
{27rduJo/cy{N/S){c/2(3) [emyc-i/v,)L _ e-2/3L/^«] 

{F}{0,t-L/vg)F,{0,t-L/vg)) 
{2nduJo/cyNL/S 



(62) 



In the denominators we have taken t = L/vg for the time over which the atoms 
radiate, and have used the fact that 2/5(l/c — l/vg)L = gL, the difference of two 
numbers that themselves are small according to our assumption that propagation 
times are small compared with the single-atom radiative decay rate, is << 1. 

The numerator in equation (|62D can be related to the expectation value q of the 
number of photons in the incident signal pulse as follows. The expectation value 
of the incident signal intensity is 

l,{0,t) = |^(F,(0,t)F,(0,t)) = V"*'/-^' (63) 

for a Gaussian pulse of duration Tp. Requiring that the energy flux j'^^dtls{z^t) 
be qhu/S = qhuo/ S implies Jq = qTioJoj {STp^) and therefore 

(Ft(0,t - L/^,)A(0,t - L/vg)) = g J^e-(*-^/^«)^/^.^ . (64) 

Thus 

Tp^TlVg \ hcb I 

A {^pLTpJ Vg 

TT Tp Vg 
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where we have used equation (|2l) [ITT 



Among the criteria given by CKK for the observation of a superluminal pulse is 
that "The probe-pulse duration [xp] must not exceed tr = ic/LUp." This criterion 
implies, from equation (p5|), that SNR{L,t) > {q/ y/7i)c/vg and therefore that it is 
possible, even for g ~ 1, to have superluminal propagation with SNR{L,t) > 1 if 
the pulse duration is short enough: Tp < TRc/vg. 

In order to relate this conclusion to ARS, we use equation (|l^ to write ( |65D as 



SNR{L, t) = e-^t-LMVr-, (66) 



We see from this expression that, if we impose the ARS condition (2), i.e., {vg/c — 
l)L/c >> Tp, then 

SNRiL, t) « S^^—e-^'-^M'l^'. , (67) 

so that, given also the condition on | A|rp discussed before equation (|5T|), the signal- 
to-noise ratio will he very small when the ARS condition for strong distinguisha- 
bility of superluminal propagation from propagation at the speed c is satisfied. 

In fact if {vg/c — l)L/c » Tp and therefore SNR{L, t) is very small for g ~ 1, 
then 

t/Tn = ^ = ^ , (68) 

Tr C Tr C Tr 

which, from (|65|), must be large. Then the SF noise must be exponentially large 



[equation (^T])]. It follows that q must be exponentially large in order to maintain a 
signal-to-noise ratio greater than unity. This is consistent with the ARS conclusion 
that "for the signal amplitude to be larger than the amplitude of the fluctuations 
at the observation time, the signal amplitude should be exponentially large" [§]. 

Our results are therefore in agreement with ARS in that, if we require the 
separation of the superluminal pulse and a twin vacuum-propagated pulse to be 
much larger than the pulse duration, the signal-to-noise ratio will be very small at 
the one- or few-photon level. On the other hand, the results are not inconsistent 
with CKK: even at the one-photon level we can achieve a signal-to-noise ratio 
greater than unity if this separation [[vg/ c — 1]L/ c) is smaller than the pulse 
duration Tp [Equation 



Physical Origin of the Noise Limiting the Observation of 
Superluminal Group Velocity 

Note that, when we set the time t in equation (^) for the short-time SF noise 
intensity equal to the "observation time" L/c, we obtain exactly the noise in- 
tensity appearing in the denominator of equation (|62D [[T^ . Thus the quantum 
noise that imposes limitations on the observation of superluminal group velocity 
is attributable to the initiation of SF. 
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Operator Ordering and Relation to ARS Approach 

Less obvious, perhaps, is the relation between the quantum noise we have consid- 
ered - whic h ste ms from the atomic dipole fluctuations characterized by equations 
(|140|) and (|141D of Appendix B - and the quantum noise of ARS, which is at- 
tributed to the quantum fluctuations of the field. 

To estabhsh the relation to the ARS approach we return to our calculation of 
the noise intensity, using now antz-normally ordered field operators instead of the 
normally ordered operators used before. Thus we consider now the expectation 
value {F{z,t)F^ {z,t)) instead of {F\z,t)F{z,t)) . In this approach the atomic 
dipole fluctuations play no explicit role, as can be seen from equation (^) and the 
fact that 

{s{z',to)s\z",to)) = Q (69) 

for excited atoms. In this case, however, the initially unoccupied modes of the field 
make a nonvanishing contribution as a consequence of non-normal ordering: 

{F{Q,t-Llv,)F\^,t-Llv,)) = E^(«'^(0)aI(0))e^(-'=)^ 

= E^-^laML + i], (70) 



which follows from (^) and ( pSD and the approximation gL << 1 upon which ( p5D 
is based. The contribution from the term that does not vanish as L — can be 
ignored, as it corresponds to vacuum quantum noise (energy ^huk per mode) that 
is present even in the absence of the amplifier. In other words, the quantum noise 
of the field in the presence of the amplifier is 

(F(0, t - z/vg)F\Q, t - z/v,))n = E -q7^^(^'^)^ ^ 2^ y du;—-g{u)L 

k 

where we have gone to the mode continuum limit, approximated uj by uJo in the 
numerator of the integrand, and used equation ( |55| ) for the gain coefficient. Per- 
forming the integration, we obtain exactly the noise term appearing in the denom- 
inator in the last line of (|6^) . But now the noise is attributable to the amplification 
of vacuum field fluctuations ||13 



Thus we can attribute the quantum noise that limits the observation of superlu- 
minal group velocity to either the quantum fluctuations of the field in the inverted 
medium, as do ARS, or to the quantum fiuctuations of the inverted atoms, as in 
our derivation of the signal-to-noise ratio. The situation here is similar to that in 
the theory of the initiation of SF, as discussed by Polder et al. [0, or, as noted by 
those authors, to the theory of spontaneous emission by a single atom M . 
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Limit of Very Small Transition Frequency 

Since the origin of noise in the optical amphfier is associated ultimately with spon- 
taneous emission, the question arises as to whether the signal-to-noise ratio might 
be increased by employing a transition having a very small transition frequency 
ojo and therefore a very large radiative lifetime. Indeed, since oc cJq, the second 
line of equation ( p3| ) suggests at first glance that SNR ^ oo in the limit Uo — >■ 0. 
However, equation ( |l6l) shows that Vg ^ c m this limit: the superluminal effect 
itself becomes weaker as the spontaneous emission rate is made smaller. 

In this connection we invoke once again the form (|66D of the signal-to-noise 
ratio. If we assume |A|rp > 1 in order that the pulse does not undergo substantial 
distortion as a consequence of strong absorption, then 

SNR{L,t)<q '\ . (72) 
{Vg -c)L/c 

In other words, the signal-to-noise ratio must be smaller than the number of 
photons in the incident pulse times a factor equal to to length of the vacuum- 
propagated pulse divided by the separation of the vacuum-propagated pulse and 
the pulse emerging from the amplifier, independent of the the atomic transition 
frequency or the radiative lifetime. At the one- or few-photon level the signal- 
to-noise ratio must therefore be less than unity under the ARS criteria for the 
observation of superluminal group velocity, regardless of the frequency or strength 
of the amplifying transition. 



5 Unitarity and Superluminal Propagation 

We now turn our attention from the specific example of the optical amplifier to 
some general features of superluminal propagation that follow generally from the 
unitary evolution of the state vector, considered here within first quantization. 

The time evolution of a wave packet can be formulated in terms of a unitary 
operator U{t) or equivalently in terms of a coordinate-space propagator G{x — 
x',t) = {x\U{t)\x'): 

\m) = u{t)\m) 

/oo 
dx'G{x-x',t)'^{x',0). (73) 
-oo 

The assumption that the propagator vanishes identically outside the light cone 
implies that 

G{x-x' > ct,t) = 0. (74) 

Given an initial wave packet centered around x = Xq < at t = 0, we assume 
that at a later time t > it will be centered around Xq + Vgt, as in the example of 
pulse propagation in an inverted medium. 
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We divide the wave packet into two parts, which we label as "superluminal" 
(S) and "luminal" (L), in the following way: 

\E''^ vanishes if the group velocity Vg < c. 

Suppose that Vg > c and that we let the wave packet propagate for a time 
T long enough that a superluminal signal can be clearly identified. That is, we 
assume that at t = T, 

^^(x,T)^0. (76) 

Now 

(^'(r)|vl>(T)) = (vl/(0)|v^(0)) (77) 

due to unitarity, and thus 

/oo roo roo 

dx'\^{x',0)\^ = dx'\^{x',T)\^ ^ dx'\^^{x',T)\\ (78) 
-oo J —oo J cT 

Physically, this means that the superluminal signal, \E'"^(x,T), is about as large, or 
contains about "as many photons," as the initial wave packet. 

We now combine the two underlying premises of causality and superluminal 
propagation as they are defined by equations (0) and ([76D. Using equation ( [73D 
for X > cT, we write 

/O rco 
dx'G{x -x',T)^{x',0) + / dx'G{x -x',T)^{x',0). (79) 
-oo JO 



The first term vanishes because, according to Eq. ([7^), the integrand differs from 
zero only ii x' > x — cT > 0. Thus 

/oo 
dx'G{x-x',T)[Q{x')'^{x',0)]. (80) 
-oo 

This formulates the notion, which is essential to the ARS argument, that for a 
causal, [i.e., equation (|7^) ], superluminal signal [equation (|76D], the wave packet is 
reconstructed from its tail [equation (80)]. 

This rather remarkable reconstruction of the signal propagated without distor- 
tion and with superluminal group velocity is especially evident in the temporal 
domain |jT5[. (See Figure 1.) 

The construction ( pOD of the superluminal wave packet from the tail of the initial 
wave packet motivated ARS to define another, truncated initial wave packet: 

$(x,0) = e(a;)^(x,0). (81) 

The two dijferent initial wave functions, ^{x,0) and \&(a;,0), give the same super- 
luminal signal: 

/oo 
dx'G{x-x',T)<!>{x',0). (82) 
-oo 
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Equation (p^ ) implies what ARS call amplification: a "small" signal propagates 
to become a "large" signal. After all, $(2;, 0) is "made from a small number of 
photons," while we have just seen that \E'(x > cT, T) has about the same number 
of photons as the non-truncated initial wave packet. We note that amplification 
in this sense is a necessary consequence of a superluminal group velocity. 

One might be tempted to write (^21) symbolically as 

|$(0)) — > l^^(r)), (83) 

where — >■ denotes time evolution under U{T). This would be incorrect: the 
truncated initial wave packet 0) is a perfectly well defined initial state, but 
it does not evolve into \E''^(x,T); part of it evolves luminally. It will be prove 
convenient to introduce "superluminal" and "luminal" parts of the truncated wave 
packet in a manner similar to the decomposition (|7^) used for the complete wave 
packet "^{x, T): 

We note that, while the superluminal part of the time-evolved truncated initial 
state is the same as the superluminal part of the time-evolved non-truncated initial 
state, the luminal parts of these signals differ: 

$^(a;,r) = ^^(x,r), (85) 
$^(x,T) ^ ^^(x,r) ^ 0. (86) 

That is, while the luminal part of the time-evolved complete wave packet approx- 
imately vanishes {^^{x,T) 0), the luminal part of the truncated wavepacket, 
$^(x, T), does not. We show below that, on the contrary, it grows exponentially 
with time. 



Momentum space: normal and unstable modes 

We are comparing the time evolution of two different initial wavepackets, \E'(x,0) 
and <I>(x, 0) where $(a;, 0) = 0(x)^(a;, 0). It is useful to define still another initial 
wave packet, 

i?(a;,0) = e(-x)^(x,0) . (87) 

Clearly, 

^(x,0) = i?(x,0) + <l>(x,0) . (88) 

After a time T, \E'(x, 0) evolves into \E'(x, T), $(a;, 0) into $(x, T), and i?(x, 0) into 
R{x,T). The time evolution is linear and 

^(x,T) = /?(x,T) + $(x,T) . (89) 
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Fourier transforming into momentum space, we define g{k), C{k), 

/oo 
dkg{k) exp[i{kx — uj^t)], 
-co 



-co 
roc 



/CO 
dk({k) exp[i{kx — Ukt)], 
-co 

/oo 
dk^{k) exp[i(/i;x — cj^t)]. 
-oo 



/^From these definitions it is straightforward to show that 

oo k — k' — IT] ^ 



/-OO 

^{k) = — dk 

27r J-oo 



27r J-oo k — k' + ir] 
where rj is an infinitesimal positive number. /^From the identity 

+ 7-^^ = -2vr.5(A: - k') 



k — k' — ir] k — k' + ir] 



it follows that 



9{k) = Cik) + . 
Equations (]90|)-(p^ can be written as well in the following way: 



■^{x,t) = / dkg{k) Mx^t), 

J— oo 

/oo 
dkg{k) 0fc(x,t), 
-oo 

POD 

R(x,t) = / dkg{k) pk{x,t), 



4Jk{x,t) = exp[i{kx -uJkt)] , 

-i f°° , exp[i{KX - uj^t)] 
(j)k{x,t) = — dn 



where 



2ti J-oo k — k — ir] 

+i exp[i{Kx -uj)] 
Pk{x,t) = — dK -— 

/TT J-oo K — k + tr] 

Ipkix, t) = 0fc(x, t) + pkix, t) . 

For t = we obtain, as required by their definitions, 

0fc(a;, 0) = Q{x)'4)k{x, 0), 
Pfc(x, 0) = Q{-x)i)k{x,Qi). 
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We now invoke the premises of causality and superluminal propagation, focus- 
ing on the ARS model involving the dispersion relation 



ojk = cVk'^ — ■ (106) 

As long as \k\ > m this dispersion relation describes normal oscillating modes. 
Unstable modes exist for \k\ < m. One might attempt to avoid the unstable 
modes altogether by choosing a g{k) that vanishes or is negligibly small for \k\ < m. 
This can be done, for example, by choosing an initial state with a Gaussian g{k), 
centered around ko and having a width Ako such that \ko ± AAj^I » m. This 
corresponds in the case of the optical amplifier to a pulse detuning large compared 
with a radiative decay rate. It turns out, however, as might be expected from 
the example of the optical amplifier, that even for such an initial wave packet 
^(x, 0) the unstable modes play an essential role in the time evolution of both the 
truncated and the residual wavepackets, $(x,t) and R{x,t), respectively. 

Consider the integrals in ( |101| ) and ( |102| ) as contour integrals in the complex k 
plane. The integrands, analytically continued into the complex k plane, each have 
a single, simple pole above or below the real k axis at k = k ±i7], and both have 
two branch points at k = ±m, which we connect with a branch cut on the line 
segment {—m,m) on the real k axis. The contour from —oo to oo should pass, 
as usual, slightly above the real k axis (at a distance smaller than rj). As shown 
below, this ensures causality according to equation ([7^). In the limit of infinite 

lim = CK, (107) 

|k|— >oo 

and on the circle at infinity, 

Kx — uj^t — ^ k(x — ct). (108) 

For X > we can therefore close the contour integral in the upper-half plane, 
whereas for x < ct we close the contour in the lower half. In both cases the 
contributions to the integral from the arcs at infinity vanish. 

Using first the residue theorem for x > ct, we see immediately that the super- 
luminal parts of the time-evolved residual and truncated wave packets satisfy 

Pfc(x > ct,t) = 0, (109) 
(f)k{x > ct,t) = exp[i{kx — Ukt)]. (110) 

These results are not surprising, as they simply reformulate equations ( |74D and 
(p2D, respectively. For x < ct, where we close the contour in the lower half-plane, 
the integral encircles the branch cut on (— m, m). After deforming the contour 
and isolating contributions from this branch cut, we use the residue theorem and 
obtain 

Pk{x < ci,t) = exp[i{kx — Ukt)] — Ijl{x,t), (HI) 
<ct,t) = ltix,t), (112) 
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where 



F,{x, t) = -l dK --M-^-cW^^-m J] 

' 2tx Jc K-k-ir] ^ ^ 

and /(^ dhi is a closed contour circling counter-clockwise the branch cut on the line 
segment (— m, m) while not circling the poles at /c ± if]. Each of the integrals, 
/^(x, t) and /fc (x, t), is dominated by a saddle point on the imaginary k axis in the 
complex K plane and exponentially grows with time. Combining terms, we obtain 

/oo 
dk g{k) Il^{x,t), (115) 
-oo 

$(x,t) = e(x-ct)^(x,t) + e(ct-x) r dk g{k) I^{x,t). (116) 

J —oo 

The integrals give exponentially growing contributions to the luminal parts of both 
the truncated and residual wavepackets. Our choice of g{k) enforces \k±iri\ > m, 
and as a result, 

dkg{k) \lt{x,t)-iax,t)] =0 (117) 



We s ee therefore that, when the residual and truncated wavepackets ( |115| ) and 
(|116|) are combined to form the complete wavepacket \E'(x,t) [Eq. (|89|) ], the expo- 
nentially growing luminal parts cancel each other. 



Discussion and implications for quantum noise 

We are studying the time evolution of three wavepackets: the complete wavepacket 
\E'(x,t), the truncated wavepacket ^{x,t), and the retarded, residual wavepacket 
R{x,t). These three wavepackets can be decomposed in two different ways. In 
equations (^)- (|9^) they were decomposed in the usual way via a Fourier transform 
at the initial time t = into normal and unstable modes. The Fourier components 
of the truncated wavepacket ({k) and the retarded wavepacket ^(fc) are related 
to the Fourier components of the complete wave packet, g{k), by equations ( p3D 
and (P^), respectively. If we choose to construct the complete wave packet from 
normal modes g{k), where \k\ 3> m, the truncated and retarded wavepackets will 
have a strong unstable- mode component in them. This was discussed by ARS, 
who pointed out that because of the unitarity of the time evolution, the unstable 
modes are accompanied by an enhancement of the quantum noise. 

In order to identify the noise in a space-time picture we employed in equations 
(p7D-( |10^ ) a less common decomposition. The difference between (p7D-( p!0^ ) and 
lies in the order of integration. Both decompositions can be derived from 




/oo 
dkg{k) exp[i{kx — Ukt)], 
-oo 
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-t 



oo 



Hx,t)^-I dq r ^p ^(P)^^P[^(g"-".^)] , (118) 
2n J-oo J-oD q — p — ii] 

R{x,t) = — dq dp "'-'^' ^^^^ . 119 

ivr J-oo J~oo q — p + trj 

Equation (pT]) describes a wavepacket made of a superposition of oscillating waves 
il)k[x,t) = exp[i{kx — ujkt)], with the momentum distribution g{k). In (|98D and 
( PUD each of these oscillating waves is replaced by a new wavefunction, (j)k{x,t) 
and pk{x,t), respectively. The weight function for the superposition forming the 
respective wavepackets remains g{k), but k has lost its meaning as a physical 
momentum. At any time, ipk = 4>k + pk- At t = 0, (pk and pk are obtained from ipk 
by truncation. At a later time t > 0, one can distinguish between two regions. In 
the superluminal region where x > ct, (pk = ipk and pk = 0. In the luminal region 
where x < ct, (pk ^ '4^k'- While ipk is everywhere a periodic wavefunction oscillating 
in space and time, (pk is in this region exponentially growing as a function of 
both t and x\ it is not oscillating in this region. In the same retarded region pk 
has a periodic oscillating component equal to ipk and an exponentially growing 
component which exactly cancels the contribution of (pk to this region. 

The three wavepackets we consider are formed by superpositions of these dif- 
ferent wavefunctions with the same weight function g{k). They evolve in time in 
the following way. In the superluminal region x > ct the oscillating wave functions 
"^k = <Pk = exp[i{kx — iJkt)], with Uk given above, combine to form a wavepacket 
moving at the group velocity Vg > c; this is the superluminal signal. In the luminal 
region x < ct the oscillating wave functions combine to cancel each other. This 
cancellation ensures the unitary time evolution of the complete wavepacket. The 
residual part of the complete wavepacket is essential for this cancellation to occur. 

Using the language of truncated wavepackets introduced by ARS, we see that 
the superluminal signal is constructed completely from the time evolution of the 
forward tail, i.e., from the time evolution of the truncated wavepacket. This trun- 
cated wavepacket evolves with time into a combination of the superluminal signal 
and an additional, exponentially increasing part in the luminal region x < ct. As 
discussed below, this additional part that grows exponentially with time can be 
expected to be accompanied by substantial quantum noise, as ARS observed using 
a different decomposition of the same truncated wavepacket. 

The new decomposition presented here therefore leads us to conclude that the 
exponentially growing noise is mostly "luminal" and will be delayed compared 
with the superluminal signal. This conclusion is consistent with the exponentially 
growing noise due to SF in the case of the optical amplifier |T^. Looking at the 



complete wavepacket, we observe that contributions from the time-evolved residual 
wavepacket will cancel in the luminal region x < ct the contributions from the time- 
evolved truncated wavepacket. However, while the signal in this region vanishes by 
the cancellation of the two exponentially growing contributions, the noise does not 



cancel - and may be very large |16|. We note that an amplification of the signal 
in the superluminal region does occur, but our new decomposition indicates that 
this amplification is mostly a result of a rather efficient constructive interference 
of oscillating wavefunctions, while the luminal parts of the time-evolved truncated 
and retarded wavepackets appear to be controlled by the unstable modes. 
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Our analysis in this section, being based on a "first-quantization" approach in 
which the wave packets are c-numbers, not operators, has not dealt explicity with 
quantum noise. However, as in the theory of the initiation of superfluorescence 0, 
the linearity of the model resulting from the approximation that there is no change 
in the atomic inversion over the time scales of interest allows a treatment of the 
operator fields as classical, fluctuating c- number fields fl^. Thus the shaded part 



of Figure la, the "tail" from which the superluminal signal evolves, becomes in 
such a treatment the truncated signal we have considered plus a fluctuating noise 
field. In the limit of a very weak incident signal pulse, the superluminal signal will 
be dominated by the noise part rather than the signal part of the tail shown in 
Figure la, and the signal-to-noise ratio will therefore be small, consistent with the 
ARS results as well as the results obtained in Section 4 for the model of an optical 
amplifier. 



6 Summary 

We have considered the effects of quantum noise on the propagation of a pulse 
with superluminal group velocity. In the case considered by CKK where an off- 
resonant, short pulse of duration Tp propagates with superluminal group velocity 
Vg in an optical amplifier, we calculated a signal-to-noise ratio SNR and found 
that, for an incident pulse consisting of a single photon, SNR << 1 under the 
condition {vg/c — 1)L >> Tp assumed by ARS for discrimination between the 
pulse propagating in the amplifier and a twin pulse propagating the same distance 
in vacuum. This result is fully consistent with the conclusions of ARS based on 
general considerations and, in particular, the reconstruction of the superluminal 
pulse from a truncated portion of the initial wave packet. However, if we impose 
the weaker condition that {vg/c — 1)L ~ Tp, then our conclusion is that SNR > 1 
is possible. However, in this case superluminal group velocity is observable in the 
arrival statistics of many photons, not per shot. 

We showed that, in the case of the optical amplifier, the quantum noise is 
attributable to the onset of superfluorescence, and could be associated either with 
the quantum fluctuations of the field, along the lines of the ARS considerations, 
or with the quantum fluctuations of the atomic dipoles. 

We then presented some general considerations based on unitarity and causality 
and introduced a new wave packet decomposition. In particular, we considered the 
"residual" wave packet in addition to the complete and truncated wave packets 
considered by ARS. This led to the conclusion that the noise is mostly luminal, 
that in the luminal region the truncated and residual signal grow exponentially but 
cancel each other as required by unitarity, but that the noise is not cancelled. For 
the case where the propagation time is large enough for the superluminal signal 
to be clearly distinguished from a twin pulse propagated at the vacuum speed of 
light, our conclusions were again consistent with ARS. 
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Appendix A 

The Heisenberg equations of motion for the Pauli operators follow from the Hamil- 
tonian (|25| ) and the commutation relations [(3"a;,o"y] = 2iaz, etc. : 

^xj = -^oCTyj , (120) 

d-yj = uJoCTxj + '^cr;,jE{zj,t) , (121) 

= -'^ayjE{zj,t) , (122) 
or, in the rotating-wave approximation, 

o-j = -iUoO-j - '^-^azjE^^\zj,t) , (123) 



2id 



[E(-)(^„t)<T,-a]i?(+)(z,,t)] . (124) 



jYiom. the formal solution of the Heisenberg equation of motion for dk{t) we 
obtain, using equation (p6|) . 



^ E^+\z,,t) + E^+\z„t) . (125) 

Here 

EW(z,t) = ^Y (^)'^^.(0)e-^-V'=^ (126) 

is the homogeneous ( "vacuum" ) solution of the Maxwell equation for the quantized 
field, while E^^\z,t) is the " source" part. Now in the mode continuum limit 
Efc^(^/27r)/rfA; = (£/27rc)/rfo., 

Ol ZTCC -^-^ Jo J -co 



Nt 
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Sc 



1=1 



f dt'ai{t')6{t' - t + [z, - Zi]/c) 



Sc 



"J ft 

/ dt'a,it')6it'-t+[z,-z,]/c) 



1=1 

Nt 



induJn . . X 27iiduj. 



-^jit) + ° ^^i^ - - Zi]/c)Oizj - Zi)e{t - [zj - Zi]/c) 



Sc Sc 

= + (127) 

Here E'^~^\zj,t) denotes the field, at the position Zj of atom j, that is produced 
by all the other atoms of the medium. 

W e now use this result, and the operator identity azjO-jit) = —crj{t), in equation 
( [123D . The result is 

aj{t) = -iuJoa,{t) - Pa,{t) - ja,jit)E^+\zj,t) , (128) 

where 

and E^+\zj,t) = E^+\zj,t) + E'^+\z.j,t). 

Similarly, using the operator identity o"j(t)(3-j(t) = ^[1 +azj{t)], we obtain from 
( P^ and QTTTl) 

&,j{t) = -2/3[l + a,,{t)] - -l-[M-\z,,t)a,it) - ajit)M+\z^,t)] . (130) 

Since the expectation value (a^) of the TLA inversion operator is p2—pi = 2p2 — 1, 
where pi and p2 are the lower- and upper-state probabilities, respectively, it follows 
that 2/5 is the radiative (spontaneous emission) decay rate: 

' 2/? = ^. (131) 



Trad She 

This is not the more familiar Einstein A coefficient for spontaneous emission, 
A = 4:\d\'^LJ^/3hc^ , because it gives the spontaneous emission rate into modes 
propagating unidirectionally with a single polarization, whereas A is the sponta- 
neous emission rate into all possible field modes in free space. In fact 1 / t^iad is the 
spontaneous emission rate implicit in much of laser theory: the coefficient X^A/Sn 
appearing in the standard expression for the gain coefficient g{i'), where A is the 
wavelength, is just I/trad times the cross-sectional area S: 

9{^) = ^(^2 - Ni)C{v) = -^{N2 - N^)SC{v) (132) 
o-n Trad 
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for (nondegenerate) upper- and lower- level population densities N2 and A^^i, respec- 
tively, an atomic lineshape function C{v), and |dp/3 = d^. 



Finally we use the definitions (^) and (|30D of F and s to obtain 



l^{t) = -z(A - tp)s^it) - '-^a,jit)Fizj, t), (133) 

a,,it) = -2/3[l + - ^[F\z„t)s,{t) - sl{t)F{z„t)] (134) 

in the rotating-wave approximation. The detuning between the TLA resonance 
frequency cUq and the central field frequency u is defined as A = — cu. Replacing 
Sj{t) and d'zjit) by s{zj,t) and az{zj,t), respectively, or s{z,t) and az{z,t) in the 
continuum limit, we obtain equations (|32D and (^). 

Appendix B 

For the initial state in which all the TLAs are in the upper state, (sj(to)) = and 
{sl{to)sj{to)) = 5ij. Then the operator 

Nt 

S = Y.rs,{to) (135) 

i=l 

satisfies 

(5) = , (136) 

Nt Nt 

{S^S) = Y.T.{^Kto)s,{Q) = NT. (137) 
i=l j=l 

In the continuum limit 



S= dzs{z,0) , (138) 

L Jo 

(5+5) = -f/ dz' dz"{s\z',Qsiz",Q) , (139) 
Jo Jo 

and we can satisfy (|136|) and ( |137|) by taking 



{s{z,to)) = 0, (140) 
{s\z',Qs{z",Q) = ^6{z'-z") . (141) 
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Figure 1: Incident (a) and transmitted (b) signals for a propagation length L 
and group velocity Vg > c. It follows from the causal connection between the two 
signals that the shaded portion of (b) is completely determined by the shaded 
portion of (a). If L(l/c — l/vg) is much larger than the pulse duration, the peak 
of the transmitted signal is reconstructed from a small tail of the incident pulse. 
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